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Abstract 

We show that the accelerated expansion of the Universe can be viewed as a crossover phe- 
nomenon where the Newton constant and the Cosmological constant are actually scaling opera- 
tors, dynamically evolving in the attraction basin of a non-Gaussian infrared fixed point, whose 
existence has been recently discussed. By linearization of the renormalized flow it is possible to 
evaluate the critical exponents, and it turns out that the approach to the fixed point is ruled 
by a marginal and a relevant direction. A smooth transition between the standard Friedmann- 
Lemaitre-Robertson-Walker (FLRW) cosmology and the observed accelerated expansion is then 
obtained, so that 0,m ~ f^A at late times. 
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I. INTRODUCTION 



One of the most important questions in modern Cosmology is to understand the origin of 
the accelerated expansion of the Universe Q|. An original approach to this problem suggests 
the possibility that the cosmological dynamics is generated by strong "renormalization group 
induced" quantum effects which would drive the (dimensionless) cosmological "constant" 
X(k) and Newton "constant" g(k) to an infrared attractive non-Gaussian fixed point 

This hypothesis has been triggered by the result of several recent investigations 
Isl 0, 3,0; 12. [kJ. which support the possibility that Quantum Einstein Gravity (QEG), 
the quantum field theory of gravity whose underlying degrees of freedom are those of the 
spacetime metric, can be defined nonperturbatively as a fundamental, "asymptotically safe" 




14 1 theory. By definition, its bare action is defined at a non-Gaussian renormalization group 



16 



.1711 a suitable fixed 

mn 



(RG) fixed point. In the framework of the effective average action 

point is known to exist in the Einstein-Hilbert truncation of theory space 5|, [7J, [l0| and in the 

ity of this approximation 
2Cj suggest that the fixed 



righer-derivative generalization j^| . Detailed analyses of the reliabi 
jyOl and the conceptually independent investigations "kIIis! 
point should indeed exist in the exact theory, implying its nonperturbative renormalizability. 

Within this framework, gravitational phenomena at a typical distance scale i = k" 1 are 
described in terms of a scale-dependent effective action T k [g^ u ] which should be thought of 
as a Wilsonian coarse-grained free-energy functional. The mass parameter k is an infrared 
cutoff in the sense that 1^ encodes the effect of all metric fluctuations with momenta larger 
than k, while those with smaller momenta are not yet "integrated out". When evaluated at 
tree level, describes all processes involving a single characteristic momentum k with all 
loop effects included. 

In Re, Q, I, has been ident.ned with the effective average action for Enc.idean q nantn m 
gravity and an exact functional RG equation for the /c-dependence of has been derived. 
Nonperturbative solutions were obtained within the "Einstein-Hilbert truncation" which 
assumes 1^ to be of the form 



T fc = (lGnGik))- 1 J d A x^g-{-R(g) + 2A(fc)}. 



(1-1) 



The RG equations yield an explicit answer for the fc-dependence of the running Newton term 
G{k) and the running cosmological term A (A;). In Refs. j^jj and j3, Q it was argued that 
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they are important for an understanding of the Planck era immediately after the big bang 



and the structure of black hole singularity. However, there are indications |24j that quantum 
Einstein gravity, because of its inherent IR divergences, is subject to strong renormalization 
effects also at very large distances. In cosmology those effects would be relevant for the 
Universe at late times, and it has been speculated that they might lead to a dynamical 
relaxation of A, thus solving the cosmological constant problem 24]. 

In Ref. j^J it was then argued that the late expansion of the Universe can be viewed as a 
renormalization group evolution near a non-Gaussian infrared fixed point (hereafter IRFP) 
where G and A become running quantities at some late time. In that work a sharp transition 
between standard FLRW cosmology and accelerated RG driven expansion was supposed to 
occur at some time, but a precise link with an underlying renormalization group trajectory 
was still missing, so that the FP was reached exactly at the transition. In spite of this 
simplifying assumption, the agreement with several high redshift observations is impressive 
[2^]. The aim of this paper is to improve the model described in Refs. P, Q| by taking 
into account the scaling evolution near the IRFP. In particular, we shall show that it is 
possible to describe the accelerated expansion of the Universe as a scaling evolution into 
an attraction basin of the IRFP. Our approach is in principle very similar to the standard 
procedure used in statistical mechanics based on the linearization of the RG equation near a 
fixed point, although we do not need to explicitly solve the RG equations. As in the case of 
QCD before the discovery of asymptotic freedom, it is possible to constrain the behaviour of 
the Gell-Mann-Low function (3(g) by means of the properties of the linearized RG flow near 
a fixed point, without explicitly solving the RG equations J^. For gravity it is possible to 
explicitly calculate the critical exponents by imposing that an allowed RG flow should be 
dynamically consistent with the Bianchi integrability condition. It turns out that the only 
possible solution predicts that the Universe approaches the IRFP along one irrelevant and 
one marginal direction. 

Section 2 studies the RG-improved Einstein equations when the energy-momentum tensor 
takes the perfect-fluid form. Section 3 studies a Bianchi I model with a multicomponent 
fluid, while concluding remarks and a critical assessment are presented in section 4. 
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II. RG-IMPROVED EINSTEIN EQUATIONS 



In the following we shall present the improved RG equation in the 3 + 1 formalism. 
Let g^y be the space-time metric with signature (— , +, +, +). A "cosmological fundamental 
observer" comoving with the cosmological fluid has 4- velocity w M = dx^/dr with = — 1, 
where r is the proper time along the fluid flow lines. The projection tensor onto the tangent 
3-space orthogonal to is h^ u = g^ v + u fl u u , with h^ v h v a = CT and h fJ, u u u = 0. We denote 
by a semicolon the standard covariant derivative and by an over-dot the differentiation with 
respect to the proper time r. The covariant derivative of w M reads as 

u^v = uj^v + + -Qhf.u - u^Uv, (2.1) 

where uj^ u = h a ^h 13 u U[ a] p] is the vorticity tensor, a^ v = h a ^hP v U(a^ — \Qh llv is the shear 
tensor, = is the expansion scalar and = u^ tV u v is the acceleration four-vector; 
square and round brackets denote anti-symmetrization and symmetrization, respectively. 
The Einstein equations read 

Rp, - \Rg»u = -Ag^ + SnGT^, (2.2) 
where A = A(x^) is the position-dependent cosmological term and G = G(x^) the position- 



dependent Newton parameter. Note that, unlike the analysis by Reuter and Weyer [27[, who 
"improve" the action functional, our right-hand side in Eq. (2.2) does not contain covariant 
derivatives of the Newton parameter nor a term describing the 4-momentum carried by G 
and A. 

The energy-momentum tensor T^ u is assumed to be covariantly conserved. For a perfect 
fluid it has the form T^ v = (p + p) u^u u + p g^ u . The conservation law T^ u - v = leads to 
mass-energy conservation 

p + Q(p + p)=0, (2.3) 

and to the equation of motion 

w" + ^— ^ = 0. (2.4) 
p + p 

The Bianchi identities require the RHS of Eq. ()2.2j) to be covariantly conserved. This 
consistency condition, together with the conservation laws f!2.3j) and ()2.4|) . provides the 
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equations for A and G, i.e. (see comments below Eq. (2.13)) 

A + 8nGp = 0, (2.5) 

hThy, - 8np hTG. v = 0, (2.6) 

by projecting along and onto the hyperplane orthogonal to These equations differ from 
the consistency condition (2.20) obtained in Ref. (27 1 from an improved action functional 
with variable G and A. 

The Raychaudhuri equation is obtained with the help of the Einstein field equations and 
of Eq. jZU), i.e. 

+ ^6 2 + 2(a 2 - uj 2 ) - + 4vrG(p + 3p) - A = 0, (2.7) 

where 2c 2 = cr^a^" and 2u; 2 = uj^uj^ . The term u M is identically vanishing for homoge- 
neous spaces. The scalar curvature of the tangent space is given by 

K, = {3) R = R + 2R^u l *u u + 2d 2 - 2u 2 - -6 2 , (2.8) 

3 

which leads, by using the field equations (|2.2j) . to the Friedmann equation 

2 

K. = 2a 2 - 2u 2 - -6 2 + 16nGp + 2A. (2.9) 
3 

In homogeneous spaces, Eqs. ()2.4|) and ()2.ti|) are identically satisfied, while the Friedmann 
equation (|2.9)L the energy conservation equation ()2.3|) and the integrability condition ()2.5|) 
constitute the evolution equations for kinematical quantities. In order to integrate them, 
the evolution equations for shear and vorticity are needed, together with the dynamical 
equations for G and A which are obtained by the RG equations. The latter are obtained 
in the Einstein-Hilbert truncation as a set of /5-functions for the dimensionless Newton 
constant and cosmological constant, g and A, 

kd k g = /3 g (g,X), kd k X = f3 x (g , A) , (2.10) 

and the link with the spacetime dynamics is provided by the so-called cut-off identification 

k = k(r,p,p,e,e,...). (2.ii) 

The dots stand for all possible physical or geometrical invariants which can act as IR regu- 
lators in the fluctuation determinant of IV 



The knowledge of the precise functional dependence in Eq. (j2.11j) would then provide a 
dynamical evolution which is consistent with the full effective action at k — 0. As it was 
explained in Ref. Q] the simple choice k oc 1/t can be justified on the ground that, if there 
are no other scales in the system when the Universe had age t, fluctuations with frequency 
greater than 1/t may not have played any role as yet, and the running must be stopped at 
k oc 1/t. On the other hand, in the radiative era in the early Universe, the energy density can 

n, 

be a better candidate to cutoff the modes [28], as is also suggested from general arguments 



based on the holographic principle [23 ]. 

At the fixed point, k in Eq. 1)2.11)1 is entirely determined by dimensional analysis, and we 
must always have a power law scaling of the type k oc 1/t, k oc p 1//4 or k oc G, for instance. 
However, for a non-Gaussian fixed point 

G=fb A = \*k 2 , (2.12) 
so that GA = g*A* and from Eqs. ()2.5)) and ()2.6j) a scaling of the type 

k = C 1/ *P 1/ \ (2-13) 

with £ = a/ A*/87rg* must always hold if no spatial variation of G and A occurs, i.e. WA^ = 
hTG. v = 0. 

It should be stressed that our analysis relies heavily on the assumption of a classical 
time evolution (2.5) in the neighbourhood of the non-trivial infrared fixed point. Such an 
assumption limits the number of relevant operators that would otherwise contribute to (2.5). 
Hence the scaling in Eq. (2.13) derives from the ad hoc assumption (2.5). It remains to be 
seen whether such an approximation is appropriate for studying the very late universe (see 
further comments at the end of section 4). 



In particular 
noticed by Ref. 



i. f)2.13j) will be satisfied in any homogeneous cosmology, as already 
I for a class of FLRW universe models, and for any matter field, so 
that it is always possible to close the system of the RG improved Einstein equations in a 
mathematically consistent way in these cases 1 . On the other hand, we shall see that the 
identification k oc 1/t is always recovered from Eq. ()2.13|) in the limit of t much larger than 
all relevant mass scales in the system. 



1 Note that there is no need to further invert the relation 1)2.13)1 with p = p(t) to determine k = k(t), since 
this map is not in general one-to-one if other dynamical and/or geometrical scales occur in the system 
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III. BIANCHI I MODEL 



Let us now consider a Bianchi I model with a multicomponent fluid. The line element 
reads as 

ds 2 = -dt 2 + a\{t)dx 2 + a 2 2 (t)dy 2 + al(t)dz 2 , (3.1) 

and in this case we have three different expansion factors. During the initial early universe 
stage, the dominant contribution in the Raychaudhuri Eq. ()2.7j) is the shear, and it is not 
a priori clear which the most physically plausible cutoff could be. On the other hand, the 
evolution equations described before are completely covariant, and Eq. (|2.13|) still holds. 
The scalar curvature JC is identically zero in Eq. (2.8); moreover, by virtue of Eqs. (2.1)- 
(2.3), (2.5), (2.7) and (2.8) the shear scalar is found to obey, even with variable G and A, 
the relation 

a 2 = S/5 6 , (3.2) 

S being a non-negative real number, and the scale S(t) being defined as S = (a^a^) 1 ^ 3 . 
The expansion scalar is then given by = 3S(t)/S(t) = 3H, H being the global Hubble 
parameter. Let us further consider a two-component perfect fluid described by p = p\ + p2 
and pi = wipi, P2 = W2P2, so that, from the conservation law (j2.3j) . we immediately get 

Mi M 2 f . 

p = pi + p 2 = - „ q - , h - „„ - 5 — , (3.3) 

M.\ and M.2 being integration constants (the factor 8n has been inserted for convenience). 
On inserting Eqs. ()2.12|) . (j2.13j) and ()3.3|) into Eq. ()2.9|) we obtain a single differential 
equation for S(t), i.e. 

Mi M 2 \V2 



S =- 3 V9^(^I+s^) +3^- (3 ' 4) 

It is interesting to discuss two particular CclSGS, 1.6. ct mixture of matter and radiation 
with no shear (the usual FLRW universe), and a stiff matter dominated universe, with shear. 
In the former case the new scale is represented by the energy density of the second fluid 
component. In the latter case the new scale present in the system is the shear. 

If we thus set wi = 1/3, W2 = 0, £ = in Eq. ()3.4j) . and the solution with S(0) = is 
given by 

1 r /3 / i 4 /3 A4i 

S^) = TT~ ^21,^ + tc) ~ t^, (3.5) 
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with t c = ^8/3y/X^g~M 3 1 /4 /M 2 . If instead M 2 = and w x 

s{t) r 3 , „. C, 



2VKg*M V v y v/3 



1 (stiff matter), we obtain 

S 2 1 1/3 



(3.6) 



for a non-vanishing shear. Note that, in this latter case, the system becomes isotropic at 
late times because 

- = S _ (3 7 ) 

6 3y/Kg*Mt + xV3 

In both cases the cutoff identification in terms of the cosmic time t is recovered in the limit 
t 3> t c or t ^> S/v^M (as can be seen by direct substitution of Eq. (|3.5jl or Eq. ()3.6|) in 
Eq. ()3.3|) and by further taking this limit), which are the additional dimensionful scales 



G = —— (g* + KiC^W) , A = & (A, + /> 2 r* 2 V 2/4 



(3.8) 



of the system, but the cutoff identification in Eq. (|2.13|) with £ = yXj&Kgl is realized 
at any time. This result suggests that this type of identification (with possibly a different 
value of £) is more powerful than the choice t oc l/k, and it can possibly be used also in the 
attraction basin of the fixed point, as we shall discuss shortly. 

In fact the IRFP hypothesis described in Refs. (3, Q assumes that the fixed point is 
non-Gaussian and attractive. We can thus describe its attraction basin by means of a two- 
dimensional subspace of irrelevant (or marginal) operators. Quite generally we can write 
g(k) = g* + hik 01 and X(k) = A* + Z^ 02 ; where 9\ > and 9 2 > are the critical exponents. 
The dimensionful G and A, by virtue of Eq. (|2.13j) . read as 

By inserting Eqs. (|3.8J) in the integrability condition (|2.5|1 . we readily find that the 
only solution valid for any value of the scaling eigenoperators hi, h 2 is obtained with 
£ = a/ (A* + h 2 )/d>TTg Jf , and the critical exponents of the Universe are given by 9\ = 2 and 
9 2 = 0. In other words, the IRFP is approached along one irrelevant and one marginal 
direction. 

From the Friedmann equation (|2.9|) . in the case of S = and M. 2 = 0, we find that the 
time evolution of the scale factor is then ruled by 

• 2 y/y/tog,M hM 

w_ 3 S^-iys + s 3w+1 ' 
The solution for a dust dominated Universe w = is thus given by 

s(t) ' 3 - 





2/3 



(3.10) 



which approaches the S oc t 4 / 3 behaviour found in Ref. j^J for large times, but reproduces 
the standard S oc t 2 ^ 3 for times much smaller than the crossover time t c ~ y/h\. For t > t c 
the Universe actually enters the attraction basin of the IRFP, which is eventually reached 
for t = oo. The deceleration parameter is given by 



Q = / ; ■=. , I 3 - 11 ) 

4(3vV8^t + 2^3h) 2 

which tends to —1/4 for t ^> t c , where Q m = Q\ = 1/2 as in the original IRFP model Q|. 

The last question we would like to answer is how to justify the domain of validity of 
our scale identification k oc p 1 / 4 in terms of a cutoff function. More precisely, we would 
like to understand whether the use of the energy density as a cutoff can be understood in 
terms of suppression of the fluctuations at a scale k oc p 1 / 4 . Let us consider, for the sake 
of simplicity, a generic self-interacting scalar theory. If we denote by V(4>) the scalar field 
potential, the energy density reads as p = <p 2 /2 + V(<f>). On the other hand, the contribution 
to the fluctuation determinant in the matter sector can easily be estimated in the proper-time 
formulation of Ref. [13], where the cutoff is realized essentially by the spectrum of the second 
functional derivative of the action. In this case the modes are cut-off at a nonvanishing mass 



scale k oc a/V"(0), the prime denoting the functional derivative with respect to 0. In order 
for the energy density to represent a meaningful cutoff, it must satisfy p oc k 4 oc V" 2 . This 
relation can be satisfied if the scalar-field evolution is dominated by the potential term. In 
fact in this case (slow-roll approximation), the kinetic term is negligible as compared to the 
potential term and one must have p oc V oc V" 2 , which is always realized for a familiar 
self-interacting scalar theory of the type V oc 4 . This result suggests that, if the dynamical 
evolution is dominated by the potential term, a scaling of the type k oc p 1 / 4 encodes the 
relevant degrees of freedom whose fluctuations of momenta greater than k are suppressed. 
Of course, this approximation becomes increasingly reliable in the infrared region, where 
only low- momentum modes are taken into account. 

IV. CONCLUDING REMARKS 

We have derived a smooth transition between standard cosmology and the IRFP modified 
cosmology described in Ref. Q which uses the scaling k oc p 1 / 4 as a cut-off identification. We 
find, for the first time in the literature, that the dynamical evolution in the neighbourhood 



of the non-Gaussian infrared fixed point can be described by means of one marginal and one 
irrelevant direction as the only possible RG trajectory which is consistent with the Bianchi 
identities. It would be therefore interesting to test this model against the Snla data, and we 
plan to investigate this issue in a separate work. 

At present, a substantial amount of literature, with detailed calculations, supports non- 
perturbative renormalizability of quantum Einstein gravity, by virtue of an ultraviolet cutoff 



at a non-Gaussian renormalization group fixed point (see Ref. 3l( and references therein). 
Hopefully, the years to come will provide at least a toy-model derivation of the improved 
Einstein equations (2.2) with running to an infrared fixed point. Until this remains an open 
problem, our original calculations are only a promising indication of the potentialities of the 
approach advocated. 

For the purpose of testing the infrared fixed-point hypothesis, the improved action func- 
tional built by some of us in Ref. [3^ might also prove useful. At a deeper level, the main 
problem one faces is as follows [^J: simple local truncations are sufficient in tre ultraviolet, 
but for k — > non-local terms should be included in the truncation ansatz for T^. Although 
we have still used a strictly local truncation, the cutoff identification k = k(x) introduces 
non-local features into the theory which, under certain conditions, are equivalent to some 
of the non-local terms in the truncated effective average action (cf. Eq. (1.1)). We refer 
the reader to the Introduction of Ref. [2^ for a discussion of the partial equivalence. 

It will be also interesting to compare our accelerated expansion of the universe with the 
one recently obtained within the framework of f(R) theories (see, for example, Ref. [33]). 
As far as the general theoretical background is concerned, we refer the reader to Ref. |3jJ 
for effective-action methods. On the experimental side, the variation of G is subject to 
experimental limits that may constrain the theory 35^ and we hope to be able to discuss 
this point as well in future work. 
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